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Fully automated calculation in fermion scattering * 

A. Lorca a 

a Deutsches Elektronen-Synchrotron, DESY, Platanenallee 6, 15738 Zeuthen, Germany 

The package alTALC has been developed for fully automated calculations of two fermion production at e + e~ 
collider and other similar reactions. We emphasize the connection and interoperability between the different 
modules required for the calculation and the external tools Diana, Form and LoopTools. Results for e + e~ — > 
ff,e + e~ are presented. 
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1. INTRODUCTION 

The proposed e + e~ international linear collider 
(ILC), is expected to open a new era of pre- 
cision physics and help understanding many of 
the fundamental questions that the discoveries of 
the forthcoming large hadron collider (LHC) will 
arise. With energies up to at least 500 GeV in 
its first phase [1] , the ILC will uniquely provide a 
level of experimental precision not reached so far 
yet in particle physics, demanding an equivalent 
understanding of the theoretical predictions and 
uncertainties which, for the two fermion produc- 
tion, should not be larger than a few pcrmill. 

In order to fulfill such theoretical duty, the cal- 
culations within the present standard model must 
include typically the next or next-to-next to lead- 
ing order perturbative corrections, depending on 
the process. The case of two fermion production 
is specially interesting since it plays a key role for 

1. luminosity monitoring (Bhabha scattering), 

2. directly measuring precise distributions and 
asymmetries, 

3. establishing limits on New Physics as con- 
tact interaction, Z', leptoquarks, or gravi- 
tons in extra dimensions, 

4. quantifying backgrounds for other pro- 
cesses. 



*Work supported in part by European's 5-th Framework 
under contract HPRN-CT-2000-00149 Physics at Collid- 
ers and by the Deutsche Forschungsgemeinschaft under 
contract SFB/TR 9-03. 



Thanks to the recent progress in computer 
technologies and software applications, particle 
physics has greatly profited from an indetermi- 
nate number of tools and packages developed for 
the specific task of automating perturbative cal- 
culations. Globally speaking, today's available 
calculational tools do not individually meet the 
next properties: 

• Performing higher order corrections beyond 
leading order. 

• Free-availability and free-based code (i.e. 
neither the package nor its software require- 
ments have to be purchased). 

• Comprehensibility, both in documentation 
and installation. 

• Flexibility, being able to deal with different 
processes, theoretical models and to allow 
some modifications. 

Next, it will be shortly described some features 
of the package alTALC showing some interesting 
immediate applications in two fermion produc- 
tion. 

2. AUTOMATION 

The tool alTALC [2] is designed with the inten- 
tion of matching the desired requirements pointed 
out above. Its goal is to perform automated per- 
turbative calculations of cross sections in high 
energy physics, restricted to 2 — > 2 fermion 
processes within the electroweak standard model 
(EWSM) and quantum electrodynamics (QED). 
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The tool integrates only frec-of-charge packages 
existing on the market. As the calculation pro- 
ceeds, in a modular fashion, it profits from an 
adaptation [3] of the Diana package [4] (based 
itself on the Qgraf [5] code) for the generation 
and analysis of Feynman graphs, the Form [6] 
language when dealing analytically with the large 
expressions in the amplitudes and, finally, the 
LoopTools [7] library (also integrating the FF 
package [8]) for the numerical calculation of the 
loop integrals. 

In order to simplify and guarantee a wide 
portability during installation, it proceeds con- 
forming to the standards of the Unix packages. 
The operations include a customization of your 
system settings via configure and the later ex- 
traction of individual packages (if needed) and 
examples by a make instruction. At this point 
the installation is finished. 

A detailed discussion of the algorithms used, 
the strategy of the calculation and internal work 
of each of the modules has been already presented 
elsewhere [2,9,10]. Here I will shortly discuss the 
working environment and the examples provided 
to the user. 

2.1. Makefile as environmental language 

Automation requires effective intercommunica- 
tion and a smart way to organize the tasks. This 
second concept has been exploited through the 
Makefile environment, providing a chain depen- 
dency between files and modules. The Make- 
file environment accomplishes the running of the 
whole process, from letting the user only the re- 
sponsability of modifying the driver file to build- 
ing the tree, loop and f ortran sections. At the 
end, the compilation of the program leading to 
the final results is achieved without further inter- 
action. 

The main directory into a process contains a 
main Makefile that organize the calculation in 
a generic footing. Under a generic unique make 
command, it calls the submakcfilcs Level .mk and 
Fortran. mk, providing them the necessary vari- 
ables to carry out the tree level, the loop level 
and the compilation of the generated code for nu- 
merical output. The Makefile will consider the 
inclusion of extra files in the main process direc- 



tory as new settings for the computation instead 
of the default ones appearing in the directories 
f orm/kitf orm3/ and diana/prg/. 

As a result of the previous call the directo- 
ries tree/, loop and f ortran are created and 
filled with the log files of the diagram genera- 
tion, symbolic amplitudes, eps drawings of the 
diagrams (for the two levels) and the complete 
Fortran programs, subroutines with the form 
factors grouped by topologies and other kinemat- 
ical functions. The elapsed time will always be 
printed in the file Makefile.log. 

The numerical executable program main . out 
is located inside f ortran and on its default run- 
ning returns a sample of differential cross sections 
for five different angles at a given energy of 500 
GeV. There are 7 printout columns: first and 
second remind the energy and scattering angle. 
Born cross section, virtual loop corrections and 
soft photon radiation appear in the columns 3 to 
5 respectively, while the sum of all them corre- 
sponds to the column number 6 (the expected re- 
sult). The last column is reserved for the squared 
of the loop amplitude serving as a rough estima- 
tor of the order of the error at next perturba- 
tive order, or main result in cases where the tree 
level is not physically possible 2 . Alternative in- 
tegrated cross sections and/or forward-backward 
asymmetries are available when the logical vari- 
ables lprintics, lprintfba become 'true'. For 
them a second group of columns show the esti- 
mated integration error following Richardson ex- 
trapolation to Romberg integration [11]. 

2.2. Examples with loop corrections 

The package contains three examples located 
in the examples/ directory which the user may 
run without further considerations for test and 
learning purposes. Two of them contain radiative 
corrections: 

• Bhabha scattering in QED (e~e + — > e~e + 
in Bhabha_QED), 

• flavour changing through neutral current 
example (e~e + — > bs in leLe.bS). 

2 As in the flavour changing neutral currents processes 
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With the purpose of giving a flavour of execu- 
tion time, the examples were run on two differ- 
ent machines, representing today's typically avail- 
able computers. The results are shown in Table 
timings and take from some seconds to minutes. 
Only when the size of the produced Fortran 
code begins to be large (i.e. more than ~200KB 
for a single subroutine to be compiled), then the 
GNU compiler gets extremely overloaded slowing 
down significantly the overall time. 

Table 1 

Typical timings for the different modules. The 
technical specifications for system 'Desk(top)': 
Intel Pentium III 853MHz cpu, 256MB RAM, 
GNU g77 version 3.3.3 compiler. For 'Lap(top)' 
are: Intel Centrino 1.5GHz cpu, 512MB RAM, 
Intel if ort version 8.1 compiler. 
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3. Selected results 

3.1. Fermion anti-fermion production 

The e~e + — ► bb and e~e + — > tt are two of the 
processes where it is expected to measure, with 
relatively high statistics, the indirect impact of 
the Higgs sector due to the large couplings in- 
volved through the top quark mass. 

Two plots are presented in Fig. ^ trying to em- 
phasize the magnitude of the corrections with in- 
creasing energy for the tt and bb production. The 
figure on top shows the expected integrated cross 
section, for the whole phase space. Continuous 
lines represent the calculation at Born level while 
the others, dotted and dash-dotted lines, account 
for the corrections with two different cuts on the 
maximum energy available for the soft-photon, 
which is related with the acollinearity of the two- 
fermions in the final state. The lower graph, scans 



the expected forward-backward asymmetry, up to 
1 TeV. Note the phenomcnological implications 
of the Z-boson threshold, inverting the asymme- 
try for the b quarks and also leading to the well- 
known resonance on the integrated cross section. 
For the top, the threshold energy for production 
is above the Z mass leading to positive asymme- 
tries. 

3.2. Massive Bhabha scattering 

Besides the phenomcnological implications of 
Bhabha scattering (e + e~ — > e + e~) in accelerator 
physics, it has also been considered here for 'cali- 
bration' in automated computations. The fact of 
comparing the massive case introduce most of the 
potential sources for discrepancies between two 
different implementations. 

In Table we present values for differen- 
tial cross sections at 500 GeV for different 
cost?. The expectations of the Born level 
and 0(a) radiative corrected values are shown 
for the package alTALC 3 , another totally in- 
dependent automated package composed out of 
FeynArts+FormCalc+LoopTools, and the 
same calculation droping out terms proportional 
to electron mass. The results are impressive since 
the agreement between both packages are far from 
the massless electron limit, and they arc limited 
by the roundoff of double-precision in the table. 
This increases from ~11 to 15 digits the previous 
agreement reported by different groups in such a 
technical comparisons [12-14], performing there- 
fore a strong check of the correctness of the cal- 
culations. 

4. Conclusions 

We have focussed on the development and im- 
mediate applicability of computer algebra and 
numerical tools to the scattering of two by two 
fermions, resulting into a major initiative for fully 
automated calculations, alTALC, involving the 
feasibility of integrating the work performed by 
others and to adapt it to some specific needs. 

The direct application of this tool led to the cal- 
culation of practically all 2 — > 2 fermion processes 

3 Quadruple precision is turned on here, setting 33 digits 
of precision. 
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of present interest to 0(a). Still these predictions 
would require further treatment to meet exper- 
imental observations, such as kincmatical cuts, 
strahlung effects or QCD corrections. 
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Figure 1. Integrated cross section (above) and Forward-backward asymmetry (below) for bb (thin lines) 
and ti (thick ones) production at high energies. The corrected results (B+l+s) are computed with two 
alternatives soft-photon maximum energy fraction, lo= 0.1 and 0.05. 



